Mon. Not. R. Astron. Soc. 000, dE] (2002) Printed 5 February 2008 (MN I*TeX style file v2.2) 



Surface density profiles of collisionless disc merger 
remnants 



Thorsten Naab^'^* and Ignacio Trujillo^''^ 

- - - ^ Insitute of Astronomy, Madingley Road, Cambridge, CB3 OH A, UK 

Universitdts-Stemwarte Miinchen, Scheinerstr. 1, D-81679 Miinc hen, Germany; 
^—•J , ^ Max-Planck-Institut fiir Astronomie, Kdnigstuhl 17, Heidelberg D-691 17, Germany 

^«-^ 1 '^School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK 



J3 



Accepted ???. Received ??? in original form ' 



ABSTRACT 

We present a detailed surface density analysis of a large sample of simulated coUision- 
. less mergers of disc galaxies with bulges (mass ratios 1:1, 2:1, 3:1, 4:1, and 6:1) and 

' without bulges (mass ratios 1:1 and 3:1). A dissipative component was not included. 

\^ \ The randomly projected remnants were fit with a single Sersic function and a Sersic 

CO ' function plus an exponential. They were classified, according to their bulge-to-total 

OO ' {B/T) ratio, either as a one-component system or as a two-component systems. In 

, general projection effects change the classification of a remnant. Only merger rem- 

|0 ■ nants of discs with bulges show properties similar to observed early-type galaxies. 

Their B/T ratios are in the range 0.2 < B/T < 1. Surprisingly, the initial mass ratio 
has a weak influence on the distributions of B/T, effective radius and Sersic index n. 
^ i' For all one-component projections (« 60% of all projections) the Sersic index distri- 

bution peaks at 3 < n < 4. However, the mass ratio is tightly linked to the properties 
of the outer exponential components which resemble pressure supported, spheroidal 
halos for 1:1 and 2:1 remnants and elongated heated discs for 6:1 remnants. We found 
^ , distinct correlations between the fitting parameters which are very similar to observed 

relations (e.g. larger bulges have lower effective surface densities). No indications for a 
correlation between the surface density profiles and other global parameters like rem- 
, nant masses, isophotal shapes or central velocity dispersions are found. The remnants 

^ • have properties similar to giant elliptical galaxies in the intermediate mass regime. A 

I binary disc merger origin for all early-type galaxies, especially the most massive ones, 

is unlikely. Observed nearby merger remnants have properties similar to the simu- 
lated remnants. They can have formed from binary disc mergers and might evolve 
into early-type galaxies within a few Gyrs. 

Key words: methods:numerical - methods: N-body simulations - galaxies: elliptical 
and lenticular, cD - galaxies: formation - galaxies: evolution - galaxies: fundamental 
parameters 
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1 INTRODUCTION 

For many decades the surface brightness profiles of elliptical 
galaxies have b een co nsidered to be well described by the 
|de Vaucoulei^ i 19481) r'^^'^-profile. However, from the early 
1990's, as the quality of photometric data improved, it be- 
came clear that the de Vaucouleurs law is not valid over the 
whole observable radius range and it became necessary to 
use a fitting function which can account for variations in the 
curvature of the liglit profile (C apac cioli 1987; Davics ct alJ 
ri988l:iCaon et al.ill^93: lBurkertllT993l') . Following these pio- 



* E-mail: naab@usm.lmu.de 



neering works, several authors have showed that the surface 
brightness profiles of elliptical galaxies are better described 
by the Sersic ( 1968) r^'^'^-profile. This is not surprising as the 
shape parameter n of the Sersic function is an additional fit- 
ting parameter. Interestingly, n is strongly correlated with 
other observed global properties of elliptical galaxies derived 
independently of the r^^" fits, su c h as: total luminosity, 
effective radius, JCaon et al.lll993l: lYoung fc Currie 1994 
ll995llJerien fc Binggelll997l:IPrugniel fc Sinlienll997^■ cen- 
tral velocity dispersion iGraham et alj|200ll : IVazdekis et alJ 
' 2004}), and also the masses of super-massive black holes 
llGraham et all 200lY The non-homology of elliptical galax- 
ies has also been shown to account for a large fraction of 
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the tilt of the Fundamental Plane jTruiillo et alJl2004f) . In 

general, more massive giant ellipticals have more concen- 
trated light profiles indicated by larger values of n which is 
most like ly a direct consequen ce of the formation process. 
Recently, Jong et al.l i2004ft have performed a two com- 
ponent decomposition (inner Sersic and outer exponential 
"disc") of a large sample of 558 early-type galaxies. They 
found strong evidence for outer exponential components in 
early-type galaxies at all luminosities. However, as they have 
pointed out, it was not possible to prove without further in- 
formation about kinematics and stellar populations whether 
this analysis has any physical significance. In particular it 
was not clear whether the exponential component repre- 
sented a classical flattened, rotationally supported thin or 
thic k disc. If the e xpone ntial components turn out to be 
real Ide Jong et all i2004) found possibly meaningful corre- 
lations between the fitting parameters such as: bulge sizes 
and disc sizes; bulge-to-disc ratio, bulge shape parameter, 
and bulge effective radius. 

It is the aim of this paper to test the significance of bi- 
nary disc mergers for the formation and evolution of early- 
type galaxies and bulges by investigating the surface density 
profiles of the merger remnants. As our progenitor discs have 
properties similar to present day disc galaxies a simulated 
merger remnant will not resemble a present day early-type 
in the sense that its stellar population is too young. The 
distribution and kinematics of the stars, however, might 
be in agreement. If the discs have merged early enough, 
z > 1, having similar properties (and there are indications 
for the existence of l arge scale discs at high redshift, see 
e.g. iLabbe et al]l2003l) their stellar populations might have 
evolved into early-type populations by now (discs might on 
average have been smaller in the past but our models are 
scale-free) and a comparison to the properties of present day 
early-type galaxies, as mentioned above, is justified. On the 
other hand with a direct comparison to observed present day 
merger remnants we can test whether nearby merger rem- 
nants originate from discs and whether they will evolve into 
early-type galaxies in the future. 

Such a sa m ple has been published by 
iRothberg fc JosephI (|2P04) who have successfully fit- 
ted Sersic-profiles to surface brightness profiles of 51 
observed merger remnants. They conclude that the 
remnants might evolve into ellipticals as they can have 
(K-band) luminosities in the range of an L* elliptical 
and photometric properties that are in principal similar 
to observed ellipticals. Most of the remnants have Sersic 
indices between 1 < n < 6 with a peak at n = 2. They find 
indications for more concentrated profiles, n > 8 for some 
remnants. However, there are still some open questions, 
e.g. they do not find a correlation of the shape parameter 
with luminosity or isophotal shape (most of the remnants 
have discy isophotes) and a kinematical investigation is 
still missing and at present it seems still difficult to argue 
that present disc merger remnants evolve into bona fide 
early-type galaxies based on observations alone. 

On the theoretical side it is now almost 30 years that it 
has been proposed that elliptical galaxies might form from 
mergers of disc galaxies. As an alternative to the even older 
monolithic collapse scenario the 'merger hypothesis' became 
more and more popular in the framework of cosmological 
models where structure forms hierarchically by mergers and 



every galaxy is expected to have experienced at least one 
(maybe more) major merger during its lifetime. The first 
numerical experiments in the early 1970s have been fol- 
lowed by more and more self-consist ent collisionless major 
merger models of disc galaxies (see iNaab fc Burkerj l2003l 
and references therein). As pure disc galaxies have much 
lower central phase-space densities than ellipticals it has 
been argued for a long time that one can not make cllig- 
tical galaxies by colli sionless me rgers of pure discs ( Ostrikeil 
ll980l: ICarlberg||l98d) . However. iHernouistl Jl993hl) "solved" 
this problem by adding sufficiently massive concentrated 
bulge components to the progenitor galaxies. The surface 
density profiles of the merger remnants have been compared 
to a de Vaucouleurs r'^^^-profile and a reasonable agree - 
ment has been found jBarne jll992l: lHernouistll99l ll993bD . 
iNaab fc Burkerj feOOsI) have found a good statistical agree- 
ment between observations of intermediate mass giant ellip- 
ticals and simulations mergers of discs with varying mass- 
ratios with respect to kinematic and photometric proper- 
ties including the isophotal shape. The most characteristic 
change with mass ratio is the amount of rotation of the 
remnants. Equal mass remnants are slowly rotating whereas 
unequal mass remnants are in general fast rota tors iBarned 
119981: iNaab et al.lll999l: iNaab fc BurkerlilioO^') . There are, 
however, disagreements in the higher order moments of the 
lino-of-sight velocity distributions for fast rotating remnants 
l^aab fc Burkcrt 2001b.l . Orbital analysis of merger rem- 
nants has shown that equal mass remnants are dominated 
by box orbits whereas u nequal mass remn ants are domi- 
nated by tube orbits fsee lJesseit et al]l2005l and references 
therein). As the stellar orbits are the backbones of elliptical 
galaxies a different orbital content might lead to different 
projected surface density profiles. 

A statistical analysis of the surface density profiles of 
merger remnants has not been performed so far. Given the 
observational evidence that the Sersic-function provides a 
better fit to observed ellipticals, we decided to investigate 
the surface density profiles of collisionless disc merger rem- 
nants in more detail and to probe whether the profiles of the 
simulated galaxies are well described by a Sersic function or 
if they can be decomposed in a Sersic part and a part with 
an exponential profile (here called "disc" for convenience) 
and whether the observed relations between the structural 
parameters can be recovered. Only recently, some groups 
have used the Sersic-function to fit simu lated merger rem- 
nant s.|GoM^e^^arcf^^^B^cei3 ll2005l) found that merger 
remnants of disc galaxies with massive bulges a re mo r e con- 
centrated than pure disc mergers, j^ournaud e t alJ (|200^ 
realized that the outer parts of remnants with mass ratios 
larger than 4:1 can show an exponential surface density pro- 
file similar to a disc galaxy. 

In this paper we present the detailed application 
of an observational ana lysis method (Truiillo ct al. 20o3; 
lAguerri fc TruiiUd l20o3) for surface brightness profiles of 
early-type galaxies to the projected surface density profiles 
of simulated galaxies. The profiles are fitted with a pure 
Sersic and a Sersic -I- exponential profile taking the full pho- 
tometric information into account (e.g. seeing, ellipticities, 
and errors). This enables us to perform a self-consistent disc 
-I- bulge decomposition. In a first step, the analysis is applied 
to a statistical set of simulated collisionless merger remnants 
of progenitor disc galaxies with bulges and mass ratios of 1:1, 
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2:1, 3:1, 4:1, and 6:1 and bulge-less progenitors with mass 
ratios of 1:1 and 3:1, for comparison. 

The paper is structured as follows: In Section |2| we re- 
view the parameters of the simulations used. The surface 
density analysis is described in detail in Section |H] with an 
exemplary application to the initial conditions and individ- 
ual simulated merger remnants shown in the appendix. In 
Section^we present the results based on the whole sample. 
The resulting parameter relations are discussed separately 
in Section]^ In Section |H| we investigate the nature of the 
outer exponential component and in Section|7]we summarise 
and discuss the results. 



2 SIMULATION PARAMETERS 

The statistical set of merger simulations of disc galaxies with 
mass ratios of 1:1, 2:1, 3:1 and 4:1 is described in detail in 
iNaab fc Burkerti l|2003») . Here we only give a brief summary. 
We used the following system of units: gravitational constant 
G=l, exponential scale length of the larger disc h = 1 and 
mass of the larger disc Md = 1. Each galaxy consists of an 
exponential disc, a spherical, non-rot ating bulge with mass 
Mi, = 1/3, a Hernquist density profile (lHernauistll990tl with 
a scale length ai, = 0.2h and a spherical pseudo-isothermal 
halo with a mas s Mh = 5.8, cut-o ff radius rc = lOh and core 
radius 7 = 1/1 The equal-mass mergers 

were calculated adopting in total 400000 particles with each 
galaxy consisting of 20000 bulge particles, 60000 disc par- 
ticles, and 120000 halo particles. The low-mass companion 
contained a fraction of l/ry (mass ratio ri — 1/1, 2/1, 3/1, 
4/1 and 6/1) the mass and the number of particles in each 
component with a disc scale length of h = \/l/rj, as ex - 
pected from the TuUy-Fisher relation JPierce fc Tullvll992D . 
The galaxies approached each other on nearly parabolic or- 
bits (i n agreement with estimate s from cosmological simula- 
tions llKhochfar fc Burk ert'"2006^') with an initial separation 
of Tsep = 30 length units and a pericenter distance of rp = 2 
length units. The initial parameters for the disc inclinations 
were selected in an u nbias ed way following the procedure 
described bv lBarne3 1^9^. We have realized 16 different 
initial orientations resulting in 16 equal-mass remnants and 
32 remna nts for every higher m ass ratio, respectively (see 
Tbl. 1 in iNaab fc BurkerdlioO^) . In addition, we have re- 
run the merger simulations with mass ratios of 1:1 and 3:1 
with bulge-less progenitor discs. This resulted in 192 simu- 
lated merger remnants that are investigated in this paper. 

All simulations have been performed using the tree-code 
VINE (Wetzstein et al., in preparation) in com bination with 
the GRAPE-5/6 special purpose hardware (iKawai et alJ 
L2OOO). The Plummer-softening for the force calculation was 
set to e = 0.05. Time integration was performed with a with 
a Leap-Frog integrator with a fixed time-s tep of dt = 0.04. 
For fu rther details on the simulations see iNaab fc Burkerli 
J2003ft . 



3 ANALYSIS OF THE SURFACE DENSITY 
PROFILE AND CLASSIFICATION 

In the classical view the surface brightness profiles of ellipti- 
cal galaxies have been assumed to follow the de Vaucouleurs 
r^'^^-profile (de Vaucouleurs 1948.) 

/(r) = J(0)e-^-«^('-/'-=''''. (1) 

/(O) is the central surface brightness, is semi-major axis 
effective radius and the factor 7.67 is chosen such that the 
effective radius encloses half the total lumin osity of the 
galaxy. In this paper we use the Sersic-function llSersicll968l) 
to fit the surface density profile of simulated merger rem- 
nants. The Sersic-function is a generalized version of the 
de Vaucouleurs-profile allowing the exponent to account for 
the varying curvature of the surface brightness profiles of ob- 
served galaxies. As we have shown in Sec. there is strong 
evidence that the curvature of the light profile is a physically 
meaningful parameter which is closely linked to other prop- 
erties of ellipticals. Throughout the paper we assume that 
the surface density of the simulated merger remnants refiect 
their surface brightness with a constant mass-to-light ratio 
along the radial profile. The Sersic-function can be written 
as 

/(r) = /(0)e-'"('■'''■■=''^^ (2) 

The three free parameters (which is one more than for the 
de Vaucouleurs-profile) are: the central surface brightness 
/(O), the effective radius r^, and the shape parameter n, 
the so called Sersic-index. The factor 6„ is a function of the 
shape parameter n, and is chosen such that the effective 
radius encloses half of the total luminosity. The exact value 
is derived from r(2n)=27(2n,b„), where r(a) and 7(a,x) are 
the gamma and the incomplete gamma functions. In case of 
n = 1 the profile is exponential and in case of n = 4 the 
profile is de Vaucouleurs. 

To systematically analyse the simulated merger rem- 
nants we created artificial images of every remnant seen from 
100 random projections by binning the central 10 length 
units into 128 x 128 pixels. For the simulation parameters 
given above this resulted in 14400 images for mergers with 
bulges and 4800 for mergers without bulges in the progen- 
itor galaxies. In an automated procedure every image was 
smoothed with a Gaussian filter of standard deviation 0.05 
length scales. The surface density profile, the isophotes and 
their best fitting ellipses were determined using a data re- 
duction package kindly provided by Ralf Bender. The sur- 
face density profile was thereafter fitted as desc r ibed b elow 
using the algorithm presented in iTruiillo et all (1200 ll) and 
lAguerri fc Truiilla (2002). As the Sersic-function fit is sen- 
sitive to the detailed shape in the innermost regions of the 
surface density profiles we followed a conservative approach 
and excluded all data points at radii rfu < 0.1 (~ 2 soften- 
ing lengths) which might be affected by the numerical force 
softening. However, using slightly smaller or larger values 
for rfit did not change the global results. 

In addition to the fit of a single Sersic-function to the 
surface density profiles we performed a decomposition (usu- 
ally called "bulge-disc decomposition") allowing for a bulge 
component described by a Sersic-function and an exponen- 
tial component like in disc galaxies llFreemanlll970tl : 

I(r) = loiacir) + iBulge(r) 
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Id{0) exp (r/Zi-D) + /s(0)e 



(3) 



The meaning of the structural parameters of the bulge com- 
ponent are the same as those described for the Sersic only 
model. /d(0) is the central surface density and /id is the 
scale length of the exponential component. The set of free 
parameters is completed with the ellipticities of the bulge e b 
and the disc eo. To avoid confusion with the observational 
notation all surface densities that result directly from the 
simulated remnants are denoted by E in the following. We 
have to note at this point that an exponential profile alone 
does not prove the existence of a disc component which in 
addition has to be flattened and supported by rotation. In 
Section |S| we will show that some remnants have exponen- 
tial components that are not discs. For convenience we keep 
the notation "disc" for the exponential component as long 
as there is not other evidence than the profile shape. 

To quantify which fraction of mass corresponds to the 
bulge and which fraction corresponds to the disc we used the 
bulge-to-total (B/T) ratio. B is the projected bulge mass 
and T = B + D is the total mass (here D is the disc mass) of 
the system. Using the parameter values obtained from the 
fit we estimated B/T according to: 



B 

T 



eB)r{2nB) 



lB{0)rlg{nB/bl"^ 



es^"ii/"n ;r(2ns)(l - Es ) + Jd (0) ( 1 - Ed) 

where r(a) is the gamma function. The frequently used 
bulge-to-disc ratio, B/D, is then given by 

B/T 



B/D = 



(5) 



1 - B/T 

In some rare cases (below 2%) the application of Eqn.^] 
directly using the parameters of the fit caused problems. For 
example, if the fitted central surface density of the exponen- 
tial scale component is below the outermost measured value 
of the remnant and the scale length ho becomes unrealisti- 
cally large (i.e. larger than the size of the remnant containing 
95% of all the particles) , we would significantly overestimate 
the disc mass and underestimate B/T. For these cases we 
assumed B/T = 1. 

The size of the BD systems was estimated by comput- 
ing the global semi-major effective radius, r^c, solving the 
equation 



B/T 
1 - B/T 



27(2n,6„(r,G/reB)'^") 
r(2n) 



A BD-decomposition of the merger remnants allowed us to 
classify the surface density profiles in two categories depend- 
ing on the bulge-to-total (B/T) ratios. Projections with 
B/T > 0.7 were considered pure bulge systems and, for sim- 
plicity reasons, the structural parameters determined from 
the Sersic-only fits were used. Projections with B/T < 0.7 
were considered two-component disc-l-bulge systems. For all 
cases we compared the reduced chi-squared of the Sersic- 
only fit and the _B_D-decomposition fit. In general the fit was 
considered sufficiently accurate if the reduced chi-squared 
value of the fit was smaller than unity. For those cases where 
both fits resulted in a reduced chi-squared value smaller 
than unity (see for example second row in Fig. IXTot we 
computed the probability 
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Figure 1. Histogram of the bulge— to-total ratios {B/T) for all 
projected 1:1, 2:1, 3:1, 4:1, and 6:1 merger remnants. The shaded 
area indicates the transformation region around B/T ^ 0.7 be- 
tween disc-|-bulge systems and pure bulge systems. The curves 
are slightly shifted for better visibility. 



(4> 



= 2j{2,r,G/hD)'l.{6) 



(7) 



that the chi-squared value is due to chance jPress et alJ 
Here 7 is the incomplete gamma function and u is 
the number of degrees of freedom. If that probability for the 
Sersic-only fit was smaller than P < 0.32 the fit was con- 
sidered as sufficiently accurate and the projected remnant 
was classified as one component with B/T — 1. Changing 
the threshold probability to P = 0.5 or P = 0.2 does not 
change the results. 

As it is the aim of this study to investigate the statisti- 
cal properties of the remnants we apply identical boundary 
conditions like weights, computation of errors etc. to all pro- 
jected images even if fine-tuning might reduce the residuals 
slightly in one or the other case. 

A potential source of concern for fitting a bulge and 
a disc to the surface brightness distribution of the galaxies 
is the degeneracy of the solution. In fact, if the fit is doing 
blindly, the algorithm would try, in some cases, to accommo- 
date the exponential component fitting the inner regions of 
the profile. To avoid this our code selects an initial guess for 
the exponential component that guaranteed that the outer 
observable points are well described by this component in 
the first step of the iteration. Even so, it could be possible 
that the final solution of the fit would provide an answer 
where the inner points of the surface brightness were de- 
scribed primarily by the exponential component. To avoid 
these cases, the algorithm checks whether the central inten- 
sity of the fitted exponential component is larger than the 
central intensity of the fitted Sersic component. If this is the 
case the fit starts again changing the initial conditions until 
a proper fit is achieved. This guarantees that the Sersic com- 
ponent is always describing the inner regions of the profile 
and that B/T is meaningful. 



4 STATISTICAL ANALYSIS 

In the appendix we show a detailed analysis of the initial 
conditions and individual examples of merger remnants. As 
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Figure 2. Distribution of properties for all projections classified as pure bulges for the projected 1:1, 2:1, 3:1, 4:1, and 6:1 merger 
remnants. The data are shown as a fraction of all projected remnants at a given mass ratio. Left: The Sersic-index n. The distributions 
show a weak dependence on mass-ratio and peak at around 3 < n < 4. Right: Effective radii which peak at 1 < r* < 1.7 for all mass 
ratios. Right: 




Table 1. Photometric properties (see Figs.0-|^ for 

mergers with bulges for different mass ratios with Icr 
deviations. 



10 



n 



Figure 3. Histogram of the Sersic-indox n for all projections. In 
contrast to Fig|2]we have plotted all results from the Sersic only 
fit independent of the classification (e.g. B/T) of the projection. 
The distributions peak at similar values but show a tail towards 
larger values of n. The observed (normalised ) distribution of local 
merger remnants iRothberg fc JoseDhll20n3) is indicated by open 
squares. 



we show there the properties of the merger remnants depend 
strongly on the respective projection angle under which the 
remnant would be observed. In this section we present the 
results of a statistical analysis of all merger remnants, with 
and without bulges in the progenitor discs, by analysing 
100 random projections of every merger remnant grouped 
according to their mass ratios. The results for mergers with 
bulges are summarized in Tab.0 



Mass ratio 


1:1 


2:1 


3:1 


4:1 


6:1 


n 


3.1 


3.56 


3.60 


3.66 


3.66 


± 


0.51 


0.59 


0.62 


0.66 


0.72 




1.19 


1.27 


1.36 


1.47 


1.60 


± 


0.16 


0.23 


0.31 


0.43 


0.51 




1.15 


1.09 


1.11 


1.16 


1.20 


± 


0.08 


0.09 


0.10 


0.08 


0.1 


riB 


1.70 


1.70 


1.66 


1.39 


1.32 


± 


0.33 


0.25 


0.56 


0.41 


0.48 




0.42 


0.38 


0.39 


0.33 


0.33 


± 


0.09 


0.07 


0.19 


0.11 


0.15 


hn 


1.32 


1.29 


1.22 


1.18 


1.19 


± 


0.24 


0.20 


0.29 


0.29 


0.18 



4.1 Progenitors with bulges 

In Fig. we show the bulge-to-total (B/T) ratios for all 
merger remnants with bulges and mass ratios of 1:1, 2:1, 
3:1, 4:1, and 6:1. The transition region from bulge+disc sys- 
tems to pure bulge systems around B/T = 0.7 is indicated 
by the shaded area. Small changes of this threshold value 
would not influence the global results as all mass ratios show 
a minimum number of projections in the transition region. 
Fig. Q shows that merging increases the bulge fraction of 
the galaxy significantly. For comparison, the B/T values of 
the initial disc do not exceed B/T = 0.20 (see appendix). 
The B/T ratios for every mass ratios cover a wide range 
from 0.2 < B/T < 1 and the distributions are remarkably 
similar. For every mass ratio we would identify pure bulge 
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Figure 4. Distribution of the properties for the projections classified as D/B systems for all mass ratios. Upper left: effective sizes; 
Upper right: Sersic-indices of the bulge components; Lower left: Effective radii of the bulge components; Lower right: Scale lengths of 
the exponential components. 



systems as well as bulge+disc systems. Using B/T = 0.7 
as the division line the percentage of pure bulge systems 
would be 66% for the 1:1, 51% for the 2:1, 54% for the 
3:1, 58% for the 4:1, and 58% for the 6:1 remnants, respec- 
tively. 1:1 remnants have the largest number of projections 
which would be classified as pure bulges. There is no trend 
for remnants with mass ratios of 2:1 and higher to be less 
dominated by bulge-only projections, e.g. 6:1 merger rem- 
nants do not have significantly more projections classified as 
disc-|-bulge systems as e.g. 2:1 remnants. However, there is 
a trend that high mass-ratio remnants classified as BD can 
have lower B/T indicating that the exponential component 
is becoming more dominant. 

In the following we focus on the properties of the projec- 
tions classified as pure bulge systems with B/T > 0.7 (the 
distributions in Figs.0-|S|are shown as fractions of all pro- 
jected remnants at a given mass ratio). The distribution of 
the best fitting Sersic-indices for the five mass ratios is shown 
in Fig. |5| The distributions are very similar for all mass ra- 
tios and peak between 3 < n < 4. The bulk of the projec- 
tions lie in the range 2.5 < n < 6 which is in general in good 
agreement with observe d values of intermediate mass giant 
elliptical galaxies (e.g. IVazdekis et al.ll2004l : iTruiillo et alJ 



120041: Ide Jong et al.l2004l and references therein) . Ellipticals 
with A^B > — 18 have on average smaller shape parameters, 
massive ellipticals with Mb < — 20 and large velocity dis- 
persions have on average shape parameters n > 4. This is 
demonstrated in Fig. 13 where show the mean values of the 



shape parameters and the effective radii in comparison to a 


sample of elliptical galaxies (Pclcticr ct al. 1990; 


Caon et alJ 


[1990. 11994I: iBineeeh & Jerien 1998: 


Gutierrez et al.. .20041). 



The numerical simulations are scale free and can therefore 
be shifted to smaller as well as larger sizes along the shaded 
region. We have assumed a scaling of 3kpc for the initial disc 
scale length which is in best agreement with observations. 
Ellipticals with larger shape (smaller) parameters which are 
in general also more massive (less massive) can not be fitted 
by our merger remnants. Interestingly, our distributions of 
the shape is very similar to the dis tributions of dark m atter 
halos in cosmological simulations jMerritt et alJl2005^ . 

For a better comparison with the sample of local merger 
remnants by .Rothb crg & Joseph (2004|) who fitted only 
Sersic profiles we show in Fig. |3 the distribution of the in- 
dices for the Sersic only fit for every projection independent 
of its B/T value. The simulated distributions peak at similar 
values and now have a prominent tail towards larger values 
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Figure 5. Distribution of effective radius r* versus shape param- 
eter n for the simulated merger remnants at different mass ratios 
(symbols with error bars from Ta b. Pa nd a sample of observed 
(dots) elliptical galaxi es jPeletie^^^p ^990t [Cagn et aljIlQQOl 
' J.I2 



ll994HBinggeh fc JerieJl998HGutierrez et al.^2004^ . We have as- 
sumed an initial disc scale length of 3kpc for a good agreement 
with the data. However, the simulations are scale free and can be 
shifted horizontally along the shaded area. 



of n. The distribution is s i milar to the local merger remnants 
from lRothberg fc Joseohl (l2004^ . However, we are missing all 
projections with n 1 and the observed distribution peaks 
at lower n (n = 2). This might reflect the fact that some pro- 
genitors had less massive bulges resulting in smaller values 
of n (see Section r4.2ll . We a lso find indices larger than n = 8 
for mass ratios 4:1 and 6:l. lRothberg fc JosephI (|200^ have 
suggested that these values result from excess light at the 
center. In our case it is those remnants that show a concen- 
trated bulge plus a very clear exponential component and 
were classified as two-component systems. 

Comparing to pre vious simulations 

iGonzalez- Garcia fc Balcellj ll2005ll have found values 
between 3 < n < 8 (they only performed Sersic-fits) for 
a small set of 1:1, 2:1 and 3:1 remnants which is in good 
agreement with our results in this mass-ratio range taking 
into account that their initial bulge was more massive (the 
bulge-to-disc ratio was 1/2). In addition, they also found 
no dependence of n on the mass-ratio. 

The distribution of the effective radii for the systems 
classified as pure bulges B /T > 0.7 is shown in Fig. |21 To 
correct for the flattening we computed the circularised effec- 
tive radius re = re\/l — e using the semi-major axis effective 
radius and the ellipticity, e, of the image. The distribu- 
tions for all mass-ratios are similar and peak at 1 < r* < 1.7 
with a tail extending to values as large as =3. Compared 
to the initial conditions the photometric sizes do, on average, 
not or only weakly increase. This result seems at first sur - 
prising, however it was already found bv lHernauis3 Jl993bl) 
that the half-mass radius of the remnant can be of a similar 
size than the initial disc. For a direct comparison we have 
computed the projected half mass radius of one o f our runs 
that is identical to "1-nrb" in iHernauistI ||l9933). In three 
different projections our values (1.80,1.2,1.15) are similar to 
theirs (1.92, 1.31, 1.18) taking into account that we have 



4-5 times more particles for the luminous components and 
7 times more particles in the halo component. 

We will now focus on all projected remnants that are 
classified as disc-f bulge systems (Fig. 2J ■ In the upper left 
we show the distribution of the global effective radii as they 
have been calculated using Eqn.|H] All remnants classified as 
DB systems have similar sizes of reO ~ 1-2, similar to the 
initial conditions. Consequently, independent of B /T and 
the initial mass-ratio all remnants have similar sizes than 
the initial discs. The distribution of the Sersic-indices and 
the effective radii of the bulge components of the BD fits 
are shown in the upper right of Fig.|3] The Sersic-indices for 
the bulge components of 1:1, 2:1 and 3:1 remnants peak at 
values of ub = 1.8. For 4:1 and 6:1 remnants the distribu- 
tions peak at slightly lower values around ub ~ 1.1 with tail 
of higher values extending to ub ~ 2.8. The effective radii 
peak around 0.2 < Ves < 0.8 for all mass-ratios. The maxi- 
mum values are about a factor four larger than the effective 
radius of the bulge component when analysed for the initial 
conditions. The scale lengths of the fitted disc components 
peak around 1.2 < hn < 1.3 which is about 30% larger than 
the initial disc scale lengths (Fig. |1J 

4.2 Bulge-less progenitors 

We have re-simulated the mergers with mass ratios 1:1 and 
3:1 using progenitor discs without bulge components to in- 
vestigate the influence of the bulge on the surface density 
proflle of the merger remnant. 

In the upper panel of Fig. |n] we show the result for 
the distribution of B/T-ratios of the projected remnants. 
The result is significantly different to the analysis of the 
remnants with bulge (see Fig. 0. Most projected remnants 
are classified as two-component systems. The distributions 
peak at B/T — 0.5 for the equal-mass remnants and B/T — 
0.35 for the 3:1 remnants, almost refiecting the mass- ratios 
of the progenitor discs. From the lower two panels in Fig.|S]it 
becomes clear that the systems classified as pure bulges have 
1.2 < n < 2.5 and therefore very shallow surface density 
profiles and, if the system is classified as a two component 
system, the Sersic index of the bulge is smaller than one. 

Interestingly, Sersic indices as smal l as n = 1 have been 
obser ved for local merger remnants jRothberg fc JosephI 
^004*) which might indicate that some progenitor discs did 
not have prominent stellar bulges (or bulges with low con- 
centration) . 

For a similar set of simulations 

iGonzalez-Garcfa fc Balcellj ||200^ fitting a Sersic only 
model for all remnants find a range of 2.4 < n < 3.2. 
There might be two origins for this discrepancy: all our 
remnants that appear more concentrated show a two 
component structure and were not classified as one compo- 
nent systems. I f wc lim it ourselves to Sersic fits only like 
[Conzalez-Garcia fc Balcelld ll2005h the indices would be in 
the range 1.2 < n < 3.1. The fact that we find smaller 
values for n might be caused by the different choice of 
initial conditions for the discs. iGonzalez-Garcia fc Balcellj 
|2005) use a truncated exponential disc which leads to 
remnant profiles that can not be fitted properly in the 
outer parts. Our remnants show very exponential profiles 
in the outer parts which arc fitted with small residuals 
leading to smaller over all values of n. Despite the differ- 
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Figure 6. Upper panel: Histogram of the bulge-to— total ratio 
(B/T) for all projected 1:1 and 3:1 bulge-less merger remnants. 
The shaded area around B/T = 0.7 divides bulge + disc sys- 
tems from pure bulge systems. Middle panel: Distribution of the 
Sersic index for all projections classified as pure bulges for the 
projected 1:1 and 3:1 merger remnants without bulges. Lower 
panel: Distribution of the Sersic index of the bulge components 
for the disc+bulge fits for the mass ratios 1:1 and 3:1. 



ences in detail the results of iGonzalez-Garcia fc Balcellj 
l|200^ and our results clearly indicate that coUisionless 
merger remnants of pure disc systems do not evolve into 
concentrated systems with surface density profiles similar 
to giant elliptical galaxies. Mergers without bulges also 
have lower phase-space densities as we show for a typical 



IX 



10^ r 



10" 





Figure 7. Phase-space densities for typical 1:1 merger remnant 
with and without a bulge in the progenitor disc. The merger with- 
out the bulge has a lower phase-space density. 



example in Fig. [7| The remnant with bulge has more 
strongly bound particles and the maximum phase space 
density is at least a factor 10 larger than for the bulge-less 
remn ant. The phase-spac e dens ities have been computed 
as in lAscasibar fc BinnevI ll2005ll using the program kindly 
provi ded by the authors. This has been suggested for a long 
time lfOstrikeil ll980l : ICarlberg|ll98^ : lHernauistllT993bh and 
has to be taken into account when interpreting simulations 
of stellar mergers of pure disc galaxies. In the following we 
only consider properties of remnants of mergers with bulges 
in the progenitor discs. 



5 GLOBAL PARAMETER RELATIONS 

In Figs. |S] and |U] we show the correlation between the disc 
scale length and the central surface density of the disc as well 
as the effective radius and the mean effective surface density 
of the bulge. In both cases there is a reasonably tight correla- 
tion in the sense that larger systems have lower surface den- 
sities. The bimodal distribution in Fig. [^arises from the fact 
that there is no continuous transition between the systems 
classified as DB systems (upper left) and systems classified 
as one-component bulges (lower right). The trend for the 
cor relations is surprisin gly similar to observed galaxies (see 
e.g. lde Jong et al.l2004l : their Figs. 6 and 7). The correlation 
between the bulge effective radii and the disc scale lengths 
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is shown by the grey dots in Fig. IIUI Systems with larger 
bulges also have larger discs. O n average hp i s a fa ctor of 
3-5 larger than VeB- In contrast, Ide .Tone et all fioO?) (open 
squares in Fig. IIUI here we have scaled the obervational data 
down by a factor 3 assuming an initial disk scale length of 
3kpc for the model disc) find that bulges in observed early- 
type galaxies are on average only a factor of 2 smaller than 
the discs. This result might reflect the fixed B /T ratio in our 
initial conditions of our sample as we find a good agreement 
with observations for small bulge sizes. The few projections 
that are found to have similar sizes in both components also 
have the largest B/T values (0.7 < B/T < 0.9). As we have 
seen in Fig.^there are only a few projections in this regime. 

In Fig. nil we show, for every projection and all mass 
ratios, B/T versus the mean effective surface density, < 
Ses >, the effective radius r^B and the Sersic index of 
the bulge component ub- The points on the vertical line, 
B/T = 1, indicate the proj ected remnants classifi ed as pure 
bulges. Observations from Ide Jong et al.l i2004ll are over- 
plotted as open diamonds. There is a weak trend for galaxies 
with a larger B/T to have a smaller bulge effective surface 
density (upper panel of Fig. IIH . In addition, galaxies with 
smaller B/T ratios also have bulges with smaller efi'ective 
radii (middle panel of Fig. lTTll . There is a good general agree- 
ment with observations but still the bulk of the simulated 
projections have small bulge sizes. There is also a correlation 
between B/T and the Sersic index of the bulge. Remnants 
with ns « 1 have small values of B/T ~ 0.3 and ns is 
increasing with increasieng B/T. A similar correlation has 
beeen found by Ide Jong et alJ (|20o3) (open diamonds) for 
their large sample of observed early-type galaxies. However, 
the observed values for ns have a larger spread and appear 
to be on average 20-30% larger for B/T < 1. The correla- 
tions for Tes and ns are much tighter for B/T < 0.7. In 
this regime t hey are also much t ighter than the observed 
correlations llde Jong et alJ l2004ft which might be due to 
the fixed bulge-to-total ratio in the initial conditions. We 
found no other significant correlation of the profile proper- 
ties with global properties of the remnants like anisotropy, 
isophotal shape, or shape of the line-of-sight velocity dis- 
tribution. In particular, there is no correlation between the 
projected central velocity dispersion of the remnants (which 
varies from a mean value of a = 0.5 for 6:1 remnants to a 
mean value of cr = 0.68 for 1:1 remnants in computational 
units) and the Sersic index of the bulge. However, this cor- 
relation has been found to b e rather tight for observed el- 
lipticals ( Graha m et alJl200ll : fVa zdckis ct al. 2004). As our 
simulations are scale free we can make the simple exper- 
iment and try to shift the velocity scale of the remnants 
to the region where the observed n shows a similar distri- 
bution peaked around n = 3.5 . Using the IVazdekis et alJ 
i2004l) data this would correspond to a small velocity range 
of 130fcm/s ^ CT ^ 170fcm/s. Galaxies with lower ct have 
smaller n, those with higher dispersion tend to have larger n. 
This would shift our more massive initial disc galaxy to the 
Milky Way regime. Taking the above considerations at face 
value we might conclude that only binary mergers of evolved 
early-type spiral galaxies with masses similar to the Milky 
Way contribute to the intermediate mass population of ellip- 
tical galaxies. Turning the argument around we can almost 
certainly exclude that the most massive ellipticals are made 
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Figure 8. Disc scale length, hjj, versus the central surface density 
of the disc, So, for all remnants with bulges and all mass ratios. 
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Figure 9. EfEcctive radius of the bulge r^s versus the mean effec- 
tive surface density of the bulge, S^a, for all projected remnants 
with bulges. If the galaxy was classified as a one component sys- 
tem we assumed r^g = r^. 



by binary mergers of disc g alaxies. This is in agreem ent with 
the kinematical analysis of lNaab fc Burker3 ll2003ll . 



6 IS THE EXPONENTIAL OUTER 
COMPONENT A REAL DISC? 

In Fig. 1121 we show the fraction of the 100 random pro- 
jections of every remnant that were classified as a two- 
component systems (the rest was classified as one component 
systems) as a function of the initial geometry of the merg- 
ers. Only for a few cases, mostly for mass ratios 1:1 and 
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Figure 10. Disc scale length, h]^, versus the effective radius of 
the bulge, reS, for all projections (grey dots) assuming the best 
fittin g disk+bulge model. Observational data for E/SO galaxies 
from Ide Jong et al.l l2nO'd) are overplotted (open squares) and 
scaled down by a factor of 3 (same scaling as in Fig. |5J corre- 
sponding to an initial model disk scale length of 3kpc. 



2:1, the remnant properties are independent of projection 
effects. Statistically we find a weak influence of the initial 
mass ratio or the initial geometry alone on the classifica- 
tion of the remnants. Typically the changes in classification 
are caused by a combination of different mass-ratios, ini- 
tial geometries and projection effects. However, the presence 
of an outer exponential is a general feature which is com- 
mon to most remnants. But what is the nature of this outer 
component? Naively one would assume that it represents a 
classical exponential disc. At least for 1:1 and 2:1 remnants 
this in terpretation is questionable as e.g. iNaab fc Burkerj 
i2003h have demonstrated that the kinematical properties 
of these remnants are very similar to slow-rotating hot stel- 
lar systems. Interestingly, if investigated in detail remnants 
of all mass ratios show similar degrees of anisotropy and the 
main kinematical difference appears t o be the amount of ro- 
tational support ( Binney. .2.005.; Burkert fc Naabll2005^ . On 
the other hand, the material in the outer regions of the rem- 
nants originates from the exponential disc component of the 
progenitor discs and the exponential profile of the remnants 
could point to this origin. 

If the exponential component of the remnants was an 
outer disc it should be significantly flattened with a high 
ellipticity (e > 0.7) if seen edge-on. Dynamically it should 
be significantly supported by rotation with the local ratio 
of the line-of-sight velocity and the velocity dispersion, v/a, 
greater than unity. Proper discs in ellipticals migh t have 
significantly larger values than that jRix et al-lllioai . how- 
ever there also exist SO galaxies which have confirmed ste llar 
disc c omponents but v/a <1 {iRlx__et al. 1992; Rubin e t alJ 
Il992l) . In case of NGC 4550 this is caused by two counter- 
rotating discs which could only have been detected by detail 
kinematical analysis. In this sense the ellipticity might be 
the stronger constraint for our first order analysis. 
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Figure 11. Upper panel: Mean effective surface density of the 
bulge versus B/T. Middle panel: Effective radius of the bulge 
versus B/T. For systems with B/T = 1 we use r^g = r^. Lower 
panel: Sersic index of the bulge versus B/T. For systems with 
B/T = 1 we use ng = n. We show all projected remnants with 
bulges in the initial condition disc s. Observational data for E/SO 
galaxies from lde Jong et al.l 1200*^) are overplotted (open squares) 
for the two lower panels and scaled down for r^B by a factor of 
3 (same scaling as in Fig.|3 which corresponds to an initial disk 
scale length of 3kpc. 
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Figure 12. Number of projections (out of 100) classified as DB 
systems for mass ratios 1:1-6:1. The numbers on the abscissa in- 
dicate individual merger geometries with initial disc orientations 
as in .Naab fc Burkcrt. I.200.1) . Table 1. 
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Figure 13. Local ratio of v/cr versus e at 1.5re for the edge-on 
projections of all merger remnants. We assume disc-like properties 
if e > 0.65 and v/a > 1. 



very high ellipticities, ei.s > 0.6, and in general f/tJi.s > 1, 
in some cases even v/ai.s > 2 which is a strong indication 
for the presence of an intact disc component that was heated 
during the merger. For comparison: the initial condition disc 
has £1.5 = 0.85 and v/ai.s ~ 5. 3:1 and 4:1 remnants have 
intermediate properties for which the two-component pro- 
jections are in good agreement with the observed ellipticals 
with exponential outer profiles an d eUiptic al-likc kinemat- 
ics as chara cterized by v/ a ~ 1 llJog fc Chi trc 200^1 . As 
proposed bv iBekkl lll99d) 6:1 and 4:1 remnants are goo d 
candidates for SO galaxies (see also iBournaud et al.ll200'3) . 
They show strong evidence for a heated disc component but 
still elliptical like surface density properties in some projec- 
tions indicating that depending on the projection the rem- 
nants would be classified as SO or elliptical galaxies. We do 
not find outer components with a flattening smaller than 
e = 0.2 which is probably due to the limited heating during 
the interaction process. Note that only the strongly inter- 
acting mergers with mass ratios of 1:1 and 2:1 lead to values 

of e < 0.45 

iBournaud et al.l ll2004l . |2005^ have analysed disc merger 
simulations including star formation and argue that ellipti- 
cals with outer exponentials can be produced by mergers 
with mass-ratios in the range 4.5:1 to 7:1. Based on our 
larger statistical sample (however, without star formation) 
we conclude that all mass ratios can lead to remnants with 
outer exponential profiles and elliptical-like kinematics. We 
can think of several possible reasons for the subtle differ- 
ences in the conclusions: either the presence of a dissipa- 
tive component and star formation chang es the properties 
of mergers with mass ratios of 1:1 to 3:1 or IBournaud et all 
feoos) have missed the exponential compon ents due to their 
simple classification criterion. In addition, IBournaud et alJ 
Il2005i) use different initial conditions and a different numer- 
ical algorithm (grid code) with a spatial resolution that is a 
factor of two lower than for our simulations. Taking all these 
differences into account the general behaviour of remnants 
with mass ratios larger than 3:1 is in good agreement. 

The orbital structure of the remnants provides addi- 
tional evidence that the structure of t he exponentia l outer 
components change with mass ratio. Ijesseit et alj ll2005l) 
have shown that at radii larger than the effective radius 
1:1 remnants have equal amounts of box (box+boxlet) and 
tube orbits. In some cases major-axis tubes are also popu- 
lated to a similar degree as minor-axis tubes supporting a 
round shape. On the other hand, 3:1 and 4:1 remnants are, 
at larger radii, clearly dominated by tube orbits. Tube or- 
bits carry the angular momentum in the systems and are 
the typical orbits in stellar discs. 



To test the properties of the exponential components of 
the remnants we measure e and v/a at a radius of 1.5reG 
were the exponential component typically becomes domi- 
nant. Typical bulges in our remnants have sizes of r^s < 
0.5reG (Fig.|4j and dominate only the center. In Fig. I13l we 
show the local values of ei.5 and v/ai.s at 1.5re for every 
remnant seen along the intermediate axis (edge-on). 1:1 and 
2:1 remnants do not show high ellipticities and are dynami- 
cally hot. In this case the exponential outer component can 
not be associated with a disc. Its pro perties are more sim i- 
lar to hot outer halos of cD galaxies l|Oraham et al.l 
Remnants with mass ratios of 6:1 on the other hand show 



7 SUMMARY & DISCUSSION 

The surface density profiles of an unbiased sample of sim- 
ulated stellar merger remnants of disc galaxies with bulges 
(mass ratios 1:1, 2:1, 3:1, 4:1, and 6:1) and without bulges 
(mass ratios 1:1 and 3:1) in the progenitor discs have been in- 
vestigated in detail. Every remnant has been analysed from 
100 random viewing angles resulting in 19200 projected rem- 
nants in total. We have fitted the profiles using a single Sersic 
function, and a Sersic function plus an exponential. Based 
on the resulting bulge-to-total ratios — systems with high 



12 Thorsten Naah & Ignacio Trujillo 



bulge-to-total ratios (B/T > 0.7) were considered as one 
component — and the quality of the fits we have classified 
the projected remnants either as one- or two-component sys- 
tems. In general the classification of every remnant changes 
with its projection and we find no strong correlation between 
the surface density properties of the remnants and the mass 
ratio of the initial discs. 

CoUisionless mergers of pure disc systems do not re- 
sult in compact remnants similar to elliptical galaxies with 
large central phase space densities. Most simulated remnants 
show a two-component structure with Sersic indices of the 
bulge ub < 1. In those cases where the profile is well fit- 
ted with a single Sersic function the profiles are inconsis- 
tent with n > 2 and therefore inconsistent with proper- 
ties of observed giant elliptical galaxies confirming e arlier 
results iOstrikeJll980t ICarlberelll986l lHernauistlll992^ . An 
additional gas component in the progenitor discs might be 
able to solve this problem. It does, however, lead to a con- 
centrated peak of stars at the c enter unless very efficient 
feedback processes are included (iMihos fc Hernauisd Il996l : 
ISprineel et all l2005t ICox et alj |2005^ . Interestingly, some 
observed near by merger remnants do sh ow exponential K- 
band profiles jRothberg fc JosephlEooil which would indi- 
cate that their progenitors did not have concentrated bulges. 
Whether those observed remnants will evolve into bona- 
fide ellipticals will depend on the amount of available gas 
and when and where it is transformed into stars. Recently, 
iRobertson et al.) (|200^ showed that gas in disc mergers can 
change the central velocity dispersionsof the remnants and 
can contribute significantly to the tilt of a the Fundamen- 
tal Plane made by merger remnants. A full analysis of the 
surface density properties is, however, still missing. 

Merger remnants of discs with bulges cover a wide range 
of bulge-to-total ratios (0.2 < B/T < 1). At every mass- 
ratio more than 50% of the projected remnants can be con- 
sidered as one-component systems with a very similar Sersic 
index distributions in the range of 2 < n < 6. For all mass 
ratios the distributions peak in the range 3 < n < 4. Sta- 
tistically the photometric sizes of the remnants are of the 
same order than the size of the initial disc, independent of 
the mass ratio. In some projections the remnants can ap- 
pear 50-100% larger, in other projections the remnants can 
even appear smaller than the initial disc. Surprisingly, the 
difference in kinematic properties for different mass-ratios 
jNaab fc Burkert- 2003) is more significant than the differ- 
ence in the properties of the surface density profiles. We do 
not find correlations between the surface density profiles of 
our remnants and global properties like velocity dispersion 
or isophotal shape. As our models are scale free we can not 
explain the observed correlations between luminosity or ve- 
locity dispersion and shape index of the whole population of 
elliptical galaxies with our models alone. However, if our ini- 
tial discs were scaled to a Milky Way-type galaxy the prop- 
erties of the remnants would be consistent with intermediate 
mass early-type galaxies in the range of —20 < Mb < — 18. 
This further supports the scenario that only part of the 
present day elliptical galaxy population can have formed 
by binary mergers of discs fiNaab fc Burkcrt 2003). More 
massive elliptical s are likely to hav e formed by a different 

Sirocess, e.g. earlv llNaab et al.l2005^ or late multiple mergers 
Weil fc Hernauis^^9^ and/or niergers of early type galax- 
ies ( e.g~ lBeI^^im2005l : iGonzalez-Garci'a fc van Albadal 



I2OO5I : iNaab et a'l]|2006l . and references therein). Of course 
the above conclusion is limited to our choice of initial condi- 
tions. Allowing for varying bulge-to-disc ratios and sizes for 
the initial discs with consistent bulge profiles might lead to 
remnants that follow the observed corr elations. A first step 
in this direction has been attempted bv lAceves fc Velazaued 
(2005). They did, however, not include bulge components 
in their initial discs which is problematic (see above) and 
a consistent study is therefore still missing. Furthermore we 
have been investigating the limiting case that the progenitor 
galaxies are gas free. It is not clear in how far the presence of 
gas and the uncertain physics of star formation and feedback 
will influence the surface density properties statistically. We 
can expect stronger changes for discs with higher gas frac- 
tions. Another study that has to be done in the future. 

More than 40% of all remnants at all mass-ratios show 
clear signatures of an outer exponential component. We have 
used additional local kinematic and photometric information 
at 1.5re// to determine the nature of this exponential com- 
ponent which is not possible by photometry alone. We found 
that the properties of the exponential component changes 
with mass ratio: For 1:1 and 2:1 remnants it resembles a 
dynamically hot spheroidal halo and for 6:1 remnants it re- 
sembles a flattened disc supported by rotation. 3:1 and 4:1 
remnants have intermediate properties. It is known that the 
violent process during 1:1 and 2:1 mergers destroys the ini- 
tial disc structure while higher mass ratio mergers are less vi- 
olent an d tend to preserve the discs of the more massive pro- 
genitor iNaab fc Burkertll2003l: iGonzalez-Garcfa fc Balcellj 
2005-5 Bo urnau^^^^T ' l2005[K This interpretation is sup- 
ported by the fact that 3:1 and 4:1 remnants are signifl- 
cantly more domin ated by tube orbit s in their outer parts 
than 1:1 remnants llJesseit et aljlioosh . We have shown that 
even collisionless 1:1 remnants, despite the violent merging, 
remember their initial state to some extend and can keep 
thei r outer exponential p article distribution. Observation- 
ally, Ide Jong et alJ ll2004ll find that most early-type galax- 
ies show an exponential outer component if a two compo- 
nent fit is performed. The physical meaning of this finding 
is, however, still uncertain. We conclude that exponential 
outer components are a generic feature of binary disc merg- 
ers. They might directly indicate that mergers involving ex- 
ponential stellar discs have played an important role dur- 
ing the formation of early-type galaxies. Of course there 
are other possibilities for the origin of outer exponential 
components. They might have formed from left-over gas 
after a disc merger jNaab fc Burkertll2001al: iBarne j l20o3 : 
ISpringel fc Hernauistl2005() . accretion of small satellites, ac- 
cretion of gas from the halo followed by star formation 
etc. As all our 1:1-4:1 remnants have kin ematic properties 
that are comparable to elliptical galaxies jNaab fc Burker3 
[200jj) we can naturally explain observations of galaxies 
with elliptical-like k inematics but expone ntial-like profiles 
lljo gfc Chitrd |20o3) . In agreement with iBournaud et alJ 
820041) we find that t he effect is strongest f or high mass- 
ratios. In contrast to iBournaud et al.l j2004l) we find that 
those systems can also form from mass ratios 1:1 - 4:1 and 
we do not see a sharp transition when changing the mass 
ratio from 4:1 to 6:1. 

We have found distinct correlations between the param- 
eters from the surface density fits: (1) the sizes of the expo- 
nential component and the bulge are correlated; (2) a sys- 
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Figure A4. Projected 'Hernauist density profile witii 

scale lengtli a (open boxes). Tlie best fitting Sersic fit (solid 
li ne) has n = 2.6 as th e projection deviates significantly from 
a ide Vaucouleurl ^^4^ profile with n = 4 (dashed line) at radii 
(r/a)V4 < 0.8. 



tem with a larger bulge-to-total ratio has on average a larger 
Sersic index, a larger bulge efltective radius and a lower effec- 
tive surface density; (3) larger bulges have a lower effective 
surface density. These results are in good agreement with 
obser vations of early-type galaxies (see e.g. |lc Joiie c t alJ 
12004) . 

If compared to nearby merger remnants the range and 
distribution of Sersic indices of the remnants ar e comparable 
to the observations jRothberg fc Josephlliooi) . In particu- 
lar, nearby remnants - similar to simulated remnants - do 
not show correlations between the surface brightness pro- 
files and other global properties like velocity dispersion or 
isophotal shape. We can therefore conclude that disc galax- 
ies are the likely progenitors of the observed remnants. They 
probably evolve into elliptical galaxies in the future. The 
time-scale for this process, however, depends on the stel- 
lar populations of the remnants and can in principle be as- 
sessed applying evolut ionary models for disc galaxies (e.g. 
iNaab fc Ostrikeill20od) to N-body models. 
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APPENDIX A: INDIVIDUAL EXAMPLES 

To demonstrate the quality and the limitations of the surface 
density analysis we show and discuss the surface density fits 
of the initial condition disc and disc-|-bulge model as well as 
of one merger remnant of every mass ratio, respectively. 



Figure A5. Variation of the fitting parameters for the initial 
condition disc+bulge system with inclination angle (0° is edge- 
on, 90° is face-on). The Sersic index of the bulge, ns, the size 
of the bulge, r^g, and B/T increase from edge-on to face-on, the 
global size, r^Q, decreases due to an apparently more dominant 
bulge component. The disc scale length does not change with 
inclination. The central surface density Iq of the disc is highest 
for the edge-on case and decreases strongly towards face-on. 



Al Initial conditions 

Fig. lAll shows the real isophotes and the best fitting ellipses, 
the surface density analysis and the residuals for the initial 
condition disc-only model as seen tilted by 5°, 45° and face- 
on. The parameters for the Sersic-only and the bulge-f disc 
fits are given in the middle panel. All three projections of the 
system are classified as one component systems with close to 
exponential density profiles. The analysis program was able 
to recover the initial condition parameters almost perfectly 
if the disc is seen face-on (last row in Fig. lAll . If the galaxy 
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Figure Al. Surface density analysis for the initial disc as used for the bulge-less mergers. Left panel: Contours of the projected surface 
density (dotted) with fitted ellipses (solid) inclined by 5° (upper plot), 45° (middle plot), and seen face-on (lower plot). Middle panel: 
Surface density profiles along the major axis of the respective projection (dots with error bars). The best fitting DB- model is shown 
by the thick solid line. The disc and bulge contribution is indicated by the thin dot-dashed and the dashed line, respectively. The thick 
dotted line shows the best fitting Sersic-only model, which is in this case close to an exponential. The shaded area indicates the radial 
range influenced by the force softening and has been excluded from the fitted data points. Right panel: Residuals for the two fits and the 
corresponding reduced chi— squared values. This is a one component system with an exponential surface density profile and for the face-on 
projection we can recover the input parameters. For n = 1 the disc scale length hp is related to the effective radius like = 1.676/i£). 
Observed in projection towards edge-on the shape parameter decreases and the size increases by 10%. 
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Figure A2. Surface density analysis as in Fig. lAll but for the initial bulge system only. Note that the projected Hernquist-bulge does 
not perfectly follow an r^/^ profile. The theoretically expected value for the effective radius, = 0.36, is similar to the fitted value. The 
SD-fit does not improve the quality of the fit in this case as the reduced chi— squared for the Sersic-only fit is already well below unity. 



is seen more edge-on the shape parameter decreases and the 
size increases by more than 10%. Note that for a Sersic- 
profile with n — 1 (exponential profile) the scale length ho 
is related to the projected effective radius like re = 1.676/io. 

In Fig. IA2I we show the analysis for the bulge com- 
ponent separately. The scale length an of the Hernquist- 
spheroid is related to its effective radius, reH, as VeH = 
afr(l -I- ^a£r)/1.33. For the model used here with a scale 
length of an = at = 0.2 this would result in VeH ~ TeB = 
0.36. This value can be recovered by the fit. We note here 



that the projected Hernquist density distribution follows a 
de Vaucouleurs like r^^'^ profile only at radii larger than half 
its scale length. If analysed over the full radial range the 
best fitting Sersic-index for the projected [HernQuist (1993) 
density distribution would be n = 2.6 (see Fig. lA4ll which is 
lower than the value of our bulge, n = 3.06, (see Fig. IA2II as 
we exclude the more flattened innermost parts from the fit. 

The analysis for the composite disc-l-bulge model is 
shown in Fig. IA3I The fitted scale length of the disc seems 
in good agreement with the initial parameters (Section 
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Figure A3. Surface density analysis as in Fig. lA3( but for the initial disc+bulge model. For the face-on projection we can nearly recover 
the input parameters. However, the size, mass fraction and curvature of the bulge is underestimated if compared to the two individual 
components (see Figs . IXTl and IX2t . The effect get stronger if the galaxy is seen more edge-on. 



independent of projection. The best agreement with the 
theoretically expected values is achieved for the face-on 
projection. The disc-to-bulge ratio D/B — 4.04 is 25% 
larger than the dynamical input value of Md/Mb = 3. The 
shape parameter of the bulge, ub = 1.94, is smaller than for 
the isolated bulge. The maximum projected size of the bulge 
is equal to the isolated bulge model. If projected nearly 
edge-on the bulge-to-total ratio drops below B/T — 0.1 
and the system appears to be disc dominated except from 
the very inner parts. The apparent size and Sersic-index 
of the bulge in the edge-on projection are significantly 
smaller. Part of the outer bulge material is counted as disc 
material to guarantee an exponential projected profile - 
which does not represent the intrinsic disc distribution for 
this projection (see Fig. lAH . The bulge can then very well 
be fitted with an exponential as only the inner part of the 
bulge is accessible. In Fig. IA5I we show how the fitting 
parameters change with disc inclination. The shape index 
of the bulge changes almost linearly from ub ~ 0.6 in the 
edge-on projection to ub ~ 2 in the face-on projection. 
The bulge is also smaller if seen edge-on. The system as 
a whole appears larger in the edge-on projection. The 
global effective radius of the initial bulge-(-disc system 
has been calculated using Eqn. |S| and is in the range of 
1.3 < TeG < 1-6. The values are lower than for the pure 
disc system (Fig. lAH due to the additional centrally 
concentrated bulge component. The disc scale length, hd, 
does not change with projection, in contrast to the disc 
only model (Fig. lAlfl . The central surface density of the 
disc, however, changes strongly with inclination, especially 



at inclinations smaller than 40° (towards edge-on). This 
effect has already been investigated with a simple model 
(see e.g. ^Graham 2001). The bulge-to-total ratio is in 
the range of 0.09 < B/T < 0.2 and is increasing from 
edge-on to face-on. The increase is strongly correlated with 
the increase of ub- These results indicate that observed 
disc-to-bulge ratios of early-type disk galaxies might be 
systematically overestimated and Sersic-indices of bulges 
might be systematically underestimated if the galaxies are 
inclined with respect to the face-on projection. This effect 
might be partly responsibl e for part of the scatte r in the 
observed B/T — n relation iAndredakis et aI]ll995^ and we 
propose to investigate projection effects of observed disc 
galaxies in more detail. It also seems difficult to reconstruct 
the underlying mass distribution of the disc and the bulge 
as the low density part of the bulge is masked by the disc. 
In this way 20% - 35% of the bulge material (depending 
on the projection) could remain undetected in observed 
disc galaxies. The above results are stable and do not 
change significantly if we vary the fitting parameters like 
weights, inner boundaries, starting values etc. However, the 
details might depend very well on our choice of the initial 
conditions. The Hernquist spheroid is a good but not the 
optimal model for real bulges and more sophisticated bulge 
models that reproduce the observed properties might be 
useful for further investigations. 
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Figure A6. Surface density analysis as in Fig. lASl but for a 1:1 merger remnant projected along its three principal axes. Based on the 
selection criteria described in Section|2lthe remnant would be classified as a pure bulge system in the first and third projections and as 
a bulge+disc system in the second projection. 
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Figure All. Same as Fig. lA6) but for a 1:1 remnant without a bulge in the progenitor disc. If observed, the remnant would be classified 
as a disc+bulge system in all projections. The system is significantly less concentrated than the 1:1 remnant with bulge fFig. IA6l . 
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A2 Individual merger remnants 

In Figs.lXBlto lXTnl we show individual examples of the sur- 
face density analysis of disc+bulge merger remnants with 
mass ratios of 1:1, 2:1, 3:1, 4:1, and 6:1, respectively. A 1:1 
and 3:1 remnant of mergers of two pure disc galaxies are 
shown in Figs. lAllI and IA12I All remnants are projected 
along the three principal axes of their moment-of-inertia ten- 
sor calculated form the 40% most tightly bound particles. 
The first projection in Fig. IA6I represents a typical system 
that is classified as a one component bulge system. The free- 
dom of adding an exponential component does not increase 
the quality of the fit. The second projection is classified as 
a disc-|-bulge system based on its bulge-to-total ratio and 
the better quality of the fit. The third projection would be 
classified as a one-component bulge system as the reduced 
chi-squared of the fit is below unity and P < 0.32 (Eqn.|7|l. 
To test for the effect of resolution we have re-simulated this 
remnant with a factor of 2 and 4 more particles and did not 
find a change in the properties exceeding the errors of the 
fitting parameters (Section IA3II . 

All three projections of the 2:1 remnant shown in Fig. 
IA7I can be fitted by a pure bulge model with sufficient ac- 
curacy. The first projection of the 3:1 remnant in Fig. IA8I 
is classified as a disc-|-bulge system whereas the other two 
projections appear to be pure bulge systems. An example 
for a system with strong indications for an outer exponen- 
tial component is the 4:1 merger remnant shown in Fig. 
IA9I The first and third projections ( "edge-on" ) clearly show 
two component systems. For the face-on projection both fits 



would result in a reduced chi-squared below one. However, 
only for the DB fit P < 0.32 and the system is classified as 
a BD system. The 6:1 remnants shown in Fig. lAlOl shows 
prominent exponential components in the "edge-on" projec- 
tions. However, if seen face-on it would still be classified 
as a single bulge. Fig. lAllI shows the surface density pro- 
file of a bulge-less 1:1 merger remnant. Due to the absence 
of the bulge the profile is significantly fiatter at small radii 
than the corresponding remnant with bulge (Fig. IA6t . In 
all the projections the remnant does not have the character- 
istics of a single component system. The inner component 
is better represented by Sersic function with a shape pa- 
rameter ub < 1, even close to a Gaussian with n = 0.5. 
A 3:1 remnant without bulge is shown in Fig. IA12I In the 
first projection the remnant is classified as a two compo- 
nent system with a very shallow inner component. The two 
remaining projections appear like one component systems, 
however, with shape indices n < 2. This remnant has very 
prolate structure at the center which is typical for remnants 
without bulges. They appear more concentrated if projected 
along this structure (first rows in Fig. IXTTl and Fig |XT2t . 

Already this first analysis clarifies that different pro- 
jections of one merger remnant can lead to different fitting 
parameters. The ranges covered in B/T and Sersic-indices 
are large and, in most cases, the classification of the merger 
remnant changes depending on the projection. It is therefore 
not possible to make any reliable statement based on indi- 
vidual projections of merger remnants and only a statistical 
analysis, as presented in Section |1] of the whole homoge- 
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Figure A13. Typical distribution of the Sersic-index (upper 
panel) and the effective radius (lower panel) for a fixed projection 
of a 1:1 merger remnant classified as a bulge for 500 independent 
realizations using the boots-trapping method. The distributions 
are narrow and the errors are well below 5%. 



nous sample of simulations will allow us to draw any global 
conclusions. 



A3 Stability of the analysis 

We used the boots- trappiuR method jHevl et al.l 11994 
iNaab fc Bu^kertll2003^ to estimate the typical errors of the 
derived fit parameters. For this analysis we created 500 re- 
sampled realizations of one projection of the original N- 
particle distribution by selecting N particles from the rem- 
nant with replacement. For every realization some particles 
might appear several times and some might not appear at 
all. Thereafter we performed a full surface density analy- 
sis for all 500 realizations and computed the mean values 
and the standard deviations of the resulting fit parameters. 
We have selected two typical merger remnants as test cases: 
a projected 1:1 remnant that was classified as a were de- 
scribed primarily by the exponential component. To avoid 
these cases, the algorithm checks whether the central in- 
tensity of the fitted exponential component is larger than 
the central intensity of the fitted Sersic bulge system and a 
projected 3:1 remnant classified as a disc -I- bulge system. 

Fig. lA13l shows the distribution of values for the Sersic- 
index and the effective radius for the 1:1 test case which 
was classified as a bulge. In both cases the distribution is 
very narrow and the error is below 5%. The distribution for 
the four fitted values — Sersic-index of the bulge, effective 
radius of the bulge, scale length of the disc and bulge-to- 
total ratio — of the 3:1 remnant classified as DB is shown in 
Fig. lA14l The errors are larger than for the 1:1 remnants but 
still only of the order of 10%. The most important result of 
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Figure A14. Boots-trapping analysis for 500 independent real- 
izations (as in Fig. lAlSl of a typical 3:1 merger remnant classified 
as a disc-|-bulge system seen from a fixed projection. The errors 
for the Sersic-index of the bulge, the effective radius of the bulge, 
the disc scale length and the bulge— to— total ratio do not exceed 
the 10% level. 



this analysis is that the variation in properties for different 
projections as shown in the previous section is real and not 
dominated by errors. The errors given here are typical and 
can be used for all following plots. 
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